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These  programs  were  originally  developed  by  the  Naval  Training  Equipment 
Center  (NTEC) , Orlando,  Florida,  for  visual  flight  simulation.  They  were  con- 
verted to  run  on  the  ETL  CDC  6400  computer,  new  input  and  output  routines  were 
developed,  and  the  radar  modeling  algorithm  was  changed  to  produce  a better 
machine  readable  rather  than  better  human  readable  scene. 

RSS  is  being  used  to  determine  the  data  base  input  requirements  and  the 
radar  modeling  algorithm  parameters  necessary  for  producing  "e.orrelatable" 
reference  scenes. 
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ABSTKACT 


PRC  Information  Sciences  Company  Report  R-1938 
User's  Manual  for  the  Reference  Scene  Software 
Steven  H.  Moritz,  October  1976  (Unclassified) 


The  U.S.  Army  Engineer  Topographic  Laboratories  (ETL)  at  Port  Belvoir  is 
presently  engaged  in  a concentrated  effort  aimed  at  developing  a methodology 
for  generating  radar  reference  scenes  from  raw  cartographic  data.  A central 
part  of  this  effort  is  the  identification  of  the  minimum  set  of  radar 
attributes  required  in  such  reference  scenes  in  order  that  they  provide 
sufficient  correlation  when  compared  with  actual  PPI  radar  images.  The 
objective  of  the  effort  described  herein  is  the  development  of  Reference 
Scene  Generation  Software  (RSS)  to  be  used  by  ETL  as  a research  tool  in 

the  development  of  the  final  reference  scene  generation  criteria. 

The  RSS  is  based  on  Digital  Radar  Landmass  Simulation  Software  (DRLMS) 
provided  by  the  Naval  Training  Equipment  Center  at  Orlando,  Florida.  The 
first  step  in  the  development  of  RSS  was  to  convert  these  programs  to  run 
in-house  on  ETL's  CDC  6400.  New  input  routines  were  written  to  permit  the 
use  of  in-housc  data  bases,  and  a new  output  routine  was  written  to  permit 
the  radar  scene  output  to  be  displayed  on  ETL's  DICOMED  plotter. 

The  second  phase  of  the  program  involved  further  modifications  of  the  pro- 
grams to  make  them  more  suitable  for  correlation  work.  Amorg  the  improve- 
ments added  were  the  capabilities  to  vary  image  resolution,  size  and 
coloring.  The  software  was  also  analysed  and  corrected  to  improve  its 
geometric  accuracy. 

Finally,  a routine  was  added  to  permit  the  incorporation  of  the  altitude 
layover  effect  into  the  reference  scenes.  This  effect  produces  a non- 
uniform  radial  translation  of  the  points  on  the  radar  image  and  its 
inclusion  is  expected  to  improve  the  correlation  obtainable  with  the 
reference  scenes. 
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PREFACE 


This  work  was  principally  performed  by  Dr.  Steven  H.  Moritz,  Planning 
Research  Corporation,  under  Contract  Number  DAAK02-75-C-0098 , Radar  Programs 
Conversion,  for  the  U.S.  Army  Engineer  Topographic  Laboratories,  Fort 
Belvoir,  Virginia,  Bruce  B.  Zimmerman,  Contracting  Officer's  Technical 
Representative. 

This  document  is  a user's  manual  for  the  Reference  Scene  Software 
(RSS)  presently  being  used  by  the  U.S.  Army  Engineer  Topographic  Labora- 
tories (ETL)  at  Fort  Belvoir,  Virginia.  It  is  not  a programmer's  manual. 
This  is  to  say,  it  is  not  intended  to  provide  sufficiently  detailed 
information  to  enable  the  reader  to  perform  software  modifications.  Rather, 
it  is  intended  to  provide  the  user  with  a basic  understanding  of  the  soft- 
ware structure  and  to  provide  all  the  information  required  to  operate  the 
system. 

To  satisfy  these  requirements,  this  document  has  been  divided  into 
four  sections.  The  first  of  these  supplies  an  overview  of  system  functions 
and  component  programs  while  the  second  examines  each  program  in  greater 
detail.  Part  three  provides  detailed  examples  of  actual  deck  structures 
for  program  operation  and  part  four  is  a compilation  of  macro- flowcharts 
for  each  major  component  of  the  system.  It  is  recommended  that  these  be 
read  in  the  order  presented.  However,  Chapter  III,  on  operating  procedure 
will  stand  on  its  own  and  may  be  used  in  this  fashion  if  an  understanding 
of  program  functions  is  not  required. 

RSS  is  based  on  DRLMS  software  supplied  by  the  Naval  Training  Equip- 
ment Center,  Orlando,  Florida.  As  the  RSS  is  presently  in  an  evolutionary 
stage,  its  final  structure  and  capabilities  will  likely  be  somewhat  more 
sophisticated  than  those  described  herein.  Supplementary  versions  for 
this  document  should  be  issued  to  keep  its  content  consistent  with  the 
changing  software  structure. 
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I.  SYSTEM  STRUCTURE 

The  radar  scene  simulation  process  employed  by  this  system  consists 
of  four  steps.  The  first  step  entails  the  preprocessing  of  the  planimetry 
data  base.  This  file  contains  the  identification  and  location  of  all  of 
the  "features"  (i.e.,  cities,  lakes,  roads,  etc.)  on  the  map  from  which 
the  scene  is  to  be  made.  As  described  below,  the  operation  of  the  RSS 
requires  that  such  features  be  described  in  a special  format,  and  it  is 
the  purpose  of  the  preprocessing  step  to  place  the  planimetry  data  in 
this  form. 

Step  2 entails  a reformating  of  the  terrain  data  base.  This  file 
contains  the  elevation  at  each  point  on  the  map  and  is  initially  obtained 
from  a modified  UNAMACE  terrain  digitization.  As  will  be  shown,  this 
format  is  incompatible  with  the  processed  planimetry  data  base  and  there- 
fore must  be  altered. 
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Step  3 involves  merging  the  planimetry  data  with  the  terrain  data 
to  produce  the  radar  scene. 

Step  4 involves  reformating  the  RSS  image  file  so  that  the  scene 
can  be  displayed  on  the  DICOMED  Scanner/Plotter  at  ETL. 

We  now  proceed  to  identify  the  various  programs  and  files  involved 
in  this  procedure.  The  descriptions  presented  in  this  section  are  cursory; 
more  detailed  discussions  will  be  presented  in  Chapter  II  of  this  document 
When  reading  the  following  paragraphs  reference  should  be  made  to  Figure 
1.1. 

1 . Software  Organization 

Proarams  RSSl  through  RSS5  may  be  viewed  as  comprising  a 
planimetry  preparation  subsystem. 

The  raw  planimetry  data  file  contains  X-Y  point  pairs  defining 
the  boundary  of  each  feature  along  with  a feature  code  identifying  the 
type  of  feature  (e.g. , rccd,  city,  lake,  etc.).  This  information  is 
in  the  ETL  flatbed  vector  digitizer  format  and  contains  pen  commands 
in  addition  to  the  map  information.  RSSl  unpacks  the  feature 
information  and  throws  away  the  pen  command  codes.  It  also  performs 
certain  initial  processing  functions  to  be  described  in  Chapter  II. 

These  functions  may  be  viewed  as  comprising  an  "enhancement" 
procedure  made  necessary  by  the  data  processing  and  accuracy  require- 
ments of  RSS. 

The  input  to  program  RSS2  consists  of  the  X-Y  pairs  describing 
each  feature.  In  RSS2  this  data  is  converted  to  an  entirely  different 
form.  Instead  of  defining  features  by  X,Y  coordinates  of  points  around 
the  perimeter,  the  feature  is  now  defined  by  horizontal  strips.  This  is 
illustrated  for  a typical  feature  in  Figure  1.2  and  Figure  1.3.  First 
the  perimeter  is  generated  in  the  form  of  short  horizontal  line  segments 
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as  illustrated  in  Figure  1.2.  This  is  done  for  both  "line"  features 
such  as  roads  and  "area"  features  such  as  lakes.  Clearly  a point 
feature  such  as  a building  can  be  displayed  in  a similar  fashion  by 
using  a single  strip  of  unit  length.  For  the  area  features,  long 
horizontal  "fill"  strips  are  generated  so  that  the  area  enclosed  by  the 
feature  perimeter  is  now  part  of  the  feature.  This  is  illustrated  in 
Figure  1.3.  Each  strip  is  defined  by  the  X,Y  coordinates  of  its  left 
end,  and  by  a delta-X  value  which  is  its  length. 

In  order  to  process  the  map  within  the  space  limitations  presented 
by  computer  memory  size  it  is  necessary  to  divide  the  map  into  small 
regions.  The  data  for  each  region  can  then  be  processed  separately. 

Program  RSS3  divides  the  planimetry  strips  generated  by  RSS2  at.  fixed 
region  boundaries.  That  is  to  say,  the  strips  are  assigned  to  their 
respective  regions  and  a strip  which  overlaps  several  regions,  say  N of 
them,  is  broken  up  into  N segments  with  each  segment  being  assigned  to 
the  region  in  which  it  lies.  In  the  process  of  accomplishing  *-his, 
the  X,Y  coordinates  of  the  left  end  of  each  strip  are  converted  to 
values  relative  to  the  lower  left  hand  corner  of  its  respective  region. 

We  have  therefore  divided  our  map  into  a series  of  smaller  maps  which 
will  eventually  be  paneled  together  to  yield  the  original. 

Program  RSS4  is  a sorting  routine  which  reorders  the  planimetry 
data  generated  by  RSS3.  The  output  of  RSS4  is  an  equivalent  data  base 
in  which  the  data  is  ordered  by  the  priorities  of  Table  1.1.  In  other 
words,  the  planimetry  strips  are  ordered  first  by  ascending  region  num- 
ber and  then  by  ascending  merge  priority.  (See  the  next  paragraph  for 
a discussion  of  the  merge  priority).  For  each  feature,  the  line  segments 
are  then  ordered  according  to  increasing  Y value,  then  according  to  in- 
creasing X value.  The  sorting  procedure  is  made  necessary  by  the  proces- 
sing requirements  of  program  RSS5. 

Program  RSS5  is  a merge  program  which  is  required  because  of  a 
subtle  aspect  of  the  data  base  generated  by  the  preceding  software. 

The  problem  is  illustrated  in  Figure  1.4.  Consider  two  features,  say 
a city  and  a lake.  Each  of  these  has  associated  with  it  a series  of 
strips  describing  its  sha^e  and  location.  Now,  if  the  lake  lies  within 
the  bounds  of  the  city,  a merge  of  the  data  is  required  since  it  is  neces- 
sary to  specify  which  feature  is  to  be  "on  top",  i.e.  the  lake  strips  must 
be  copied  over  the  city  strips  or  the  lake  will  not  appear  on  the  final 
picture.  The  output  of  RSS5 , then,  is  a data  base  in  which  the  contents 
of  each  region  are  described  by  horizontal  strips,  none  of  which  overlap. 
Strips  from  features  with  low  merge  priorities  (such  as  cities)  are  entered 
first  and  then  written  over  where  necessary  by  strips  from  features  with 
higher  merge  priorities  (such  as  lakes) . The  merge  priorities  are  assigned 
by  RSS3  on  the  basis  of  feature  code. 

Program  RSS5  marks  the  end  of  the  so-called  planimetry  preparation 
subsystem. 
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Program  RSS6  is  tfclg  terrain  preparation  program.  As  indicated 
above,  the  planimetry  data  is  broken  into  regions  in  order  to  simplify 
processing.  The  terrain  data  must  be  similarly  subdivided,  and  this  is 
done  in  RSS6. 

Present  operation  of  the  simulator  calls  for  the  generation  of  four 
(4)  scenes  for  each  target.  Each  of  these  scenes  depicts  the  radar  return 
of  the  terrain  when  viewed  from  a different  altitude.  In  order  to  obtain 
such  a complete  set  of  scenes  it  is  necessary  to  run  RSS7  and  RSS8  once 
at  the  highest  altitude.  Programs  RSS9  and  RSS10  must  be  run  once  for 
each  altitude.  The  generation  of  a set  of  four  scenes  for  a given  target 
therefore  requires  one  run  of  programs  RSS1  through  RSS8  and  four  runs 
of  RSS9  and  RSS10. 

Program  RSS7  merges  the  planimetry  and  terrain  files  for  those 
regions  lying  within  the  radar  ground  range  of  the  target.  That  is  to 
say,  th=*  output  from  RSS7  consists  of  one  record  for  each  region  with- 
in the  radar  range  of  the  target.  This  record  consists  of  an  appropriate 
header,  followed  by  the  planimetry  content  of  each  point  in  the  region 
followed  by  the  elevation  at  each  point  within  the  region.  A discussion 
of  the  size  of  the  regions  and  the  resolution  cells  within  each  region 
will  be  given  in  Chapter  II. 

Program  RSS8  uses  the  cartesian  data  base  output  from  RSS7  to 
construct  a series  of  radial  scan  lines  required  to  simulate  a Plan- 
Position-Indicator  (PPI)  radar  scene  as  illustrated  in  Figure  1.5  This 
operation  consists  of  a cartesian- to-polar  coordinate  transformation. 

The  output  from  RSS8  therefore  consists  of  N scan  lines  where  360  /N 
is  the  angular  increment  between  scan  lines.  The  single  record  for 
each  scan  line  contains  the  planimetry  and  elevation  data  for  all  points 
lying  along  the  scan  line,  beginning  at  the  target  location  and  ending 
20  nautical  miles  out.  The  planimetry  data  from  RSS8  consists  of  in- 
trinsic radar  strength-of-returns  that  are  assigned  to  each  feature  on 
the  basis  of  feature  code.  Ground  points  not  assigned  to  any  special 
feature  type  are  given  a predetermined  background  return  strength. 

Program  RSS9  applies  the  radar  effects  to  the  radial  scan  line  data 
from  RS38.  At  the  present  time  only  three  effects  are  considered.  One 
is  a Lambert's  Law  effect  which  determines  the  percentage  of  the  incident 
radar  signal  reflected  back  to  the  source.  The  second  is  a shadow  effect 
which  takes  into  account  the  fact  that  certain  areas  on  the  ground  may  be 
blocked  from  view  by  tall  objects  (i.e.  mountain  peaks)  lying  in  the  line- 
of-sight  from  the  radar  location  to  the  area  in  question.  The  third  is 
the  altitude  layover  effect  which  accounts  for  the  fact  that  ground  points 
with  elevations  different  from  that  of  the  target  (i.e.  the  point  directly 
below  the  radar)  will  appear  shifted  from  their  true  position. 

These  effects  will  be  discussed  in  more  detail  in  Chapter  II.  RSS9 
also  contains  a subroutine  used  to  scale  the  final  image  so  that  scenes 
generated  from  different  altitudes  will  all  oe  the  same  size  when  displayed 
on  the  DICOMED  plotter. 
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ORDERING  1’RLORITY 

PARAMETER 

1 

REGION  NUMBER 

2 

MERGE  PRIORITY 

jf 

Y COORDINATE 

4 

X COORDINATE 

TABLE  1.1  - ORDERING  PRIORITY 

FOR  RLMS4  SORT  ROUTINE 

REGION  N 

^!Ti  irn^ 


CITY 


REGION  N 


REGION  N 


^V\ 

LAKE 


LAKE  COMPOSITE 

FIGURE  1.4  - MERGE  PROCESS  ILLUSTRATED 


FIGURE  1.5  - 


RADIAL  SCAN  LINES.  THE  RADAR  TARGET  IS 
LOCATED  AT  THE  APEX 
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At  this  point  the  radar  scene  is  completed.  The  SORT  program  and 
RSS10  are  now  used  to  convert  the  scene  information  back  to  cartesian 
form  and  to  place  the  information  on  tape  in  a raster  format  compatible 
with  the  DICOMED  plotter. 

2.  File  Structure 


The  file  structure  for  RSS  is  as  indicated  in  Figure  1.1. 

The  numbers  associated  with  each  file  in  the  diagram  are  the  file  I.D.'s 
used  in  the  programs  and  this  information  will  be  of  use  in  Chapter  III 
when  we  discuss  operating  procedures. 

As  is  evident,  all  intermediate  storage  is  done  on  disc.  Inputs  to 
the  simulator  consist  of  the  planimetry  data  tape  to  RSS1,  the  terrain 
data  tape  to  RSS6,  and  card  input  to  RSS1,  RSS2,  RSS7,  and  RSS9.  These 
card  inputs  will  be  discussed  in  Chapter  II. 


II.  SOFTWARE  DESCRIPTION 
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1.  Program  RSS1 

1.1  Inputs 

This  program  requires  both  tape  and  card  input.  The  tape 
input  consists  of  a digitizer  tape  containing  the  (X.Y)  coordinates 
of  the  points  defining  the  location  of  the  planimetry  features  on  the 
map.  The  format  of  this  tape  is  presented  in  Figure  II. 1. 

Each  small  box  represents  a 6-bit  BCD  character,  and  each  physical 
record  contains  1800  such  characters  (or  180  CDC  words).  The  beginning 
of  the  information  for  each  feature  is  marked  by  a "(".  Following  this 
are  twelve  special  characters  which  are  of  no  use  in  RSS  and  are  dis- 
carded. The  next  five  characters  contain  the  feature  code.  A 
tabulation  of  the  possible  feature  codes  is  given  in  Table  II. 1. 

Following  the  feature  code  is  a pen-down  command  which  signals  the 
beginning  of  the  (X,Y)  data. 

The  (X , Y ) data  for  the  feature  begins  in  the  twentieth  character 
of  the  feature  record.  Each  coordinate  is  in  an  F7.3  format.  This  means 
that  each  digit  (as  well  as  the  decimal  point)  is  represented  as  a 
separate  character,  with  the  maximum  coordinate  value  being  999.999. 

All  coordinates  are  given  in  inches  relative  to  the  table  origin.  The 
last  Y-value  for  the  feature  is  followed  by  a "U"  which  is  the  pen-up 
command . 

A feature  record  may  begin  anywhere  in  a block  of  data,  and  may 
require  more  than  one  block.  Similarly,  a single  block  may  contain  the 
information  for  several  features  of  small  size.  A feature  record  may 
be  broken  at  any  point  except  in  the  middle  of  the  string  of  characters 
representing  the  feature  code  or  a coordinate.  A "B"  is  used  to  signal 
the  end-of -information  for  a particular  block. 

The  card  input  to  RSS1  consists  of  a single  card  containing  the 
X and  Y offsets.  These  numbers  are  subtracted  from  every  (X , Y ) pair 
for  every  feature  and  are  required  to  bring  the  origin  of  the  planimetry 
data  into  coincidence  with  the  origin  used  for  the  terrain  data.  The 
information  is  entered  in  a 2F10.3  format. 

1. 2 Description  of  Processing 

For  each  feature,  RSS1  begins  by  pulling-off  and  recon- 
strucing  the  feature  code  and  (X,Y)  data  pairs.  All  other  information 
on  the  data  tape  is  discarded.  The  X and  Y offsets  are  substracted  from 
the  (X,Y)  pairs,  and  any  point  having  a negative  coordinate  is  discarded 
since  it  lies  off  of  the  adjusted  map  sheet.  As  mentioned  earlier,  this 
procedure  is  necessary  to  bring  the  coordinate  origin  of  the  planimetry 
map  into  coincidence  with  that  of  the  terrain  map.  Failure  to  do  this 
will  cause  a systematic  error  in  the  placement  of  the  planimetry  data. 
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FEATURE  tyre 

FEATURE  CODE 

RADAR  GREY  SHADES 

20101  THRU  20320 

LARGE  RIVERS  ( WA  TER  ON  LEFT) 

10110 

LARGE  RIVERS  (WATER  ON  RIGHT) 

1012  3 

OAKS 

10133 

MARSHES  AND  SWAMPS 

10140 

LAKES 

1 0 iso 

ISLANDS 

10160 

RIVERS  AND  STREAMS 

10170 

RAILROAD  YARDS 

10210 

RAILROADS 

10220 

TOWNS  AND  SU3UR3S 

10310 

MEDIUM  CITIES  AND  COMMERCIAL  AREAS 

10320 

BIG  CITIES  AND  INDUSTRIAL  AREAS 

10330 

LARGE  ISOLATED  3UIL0INGS 

10340 

INTERSTATE  HIGHWAYS  AND  TURNPIKES 

1 OLIO 

MAJOR  ROADS 

104  20 

SECONDARY  ROADS 

10430 

UNPAVED  ROADS  AND  TRAILS 

10440 

AIRPORT 

10450 

POWER  LINE  TOWERS  (WITH  CABLES) 

10510 

DRIVE-IN  MOVIES 

10520 

FIRE  OR  RADIO  TOWERS 

10530 

CEMETERIES 

10540 

POL  AREA 

1 0550 

HARDWOOD  FOREST 

10610 

EVERGREEN  rOREST 

10620 

MEADOWS  AND  GRASSY  FIELOS 

10630 

DRY  ROCKY  AREAS 

10640 

SANO  AND  SAGEBRUSH  AREAS 

10650 

SNOW  COVERED  AREAS 

10660 

ORY  RIVEROEJS,  CANALS,  AND  STORM  DRAINS 

1 0670 

DRY  LAKE  BEDS  AND  GULCHES 

10680 

tajujj  ij.1  - r()xio:si'ONni:N(:i:  ictiikn  feature  tyres 

* 1 • AND*  FEATURE*  (.‘CN)KS 
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The  data  for  areal  features,  such  as  lakes  or  cities,  is  then 
checked  to  see  that  a closed  polygon  is  formed;  and  closure  is  performed 
if  required.  Such  a closure  of  all  polygons  is  required  by  RSS2.  Poly- 
gons lying  partially  off  the  map  sheet  are  closed  along  the  sheet  boundary. 

In  order  to  simplify  later  processing,  the  feature  codes  of  Table 
III.l  are  mapped  into  the  set  of  integers  (1,60).  This  permits  easy 
formulation  of  translation  tables  which  match  feature  codes  with  feature 
attributes  such  as  merge  priority  and  intrinsic  radar  return  intensity. 

Finally,  RSS1  does  special  processing  on  certain  features.  If  a 
feature  is  a power  line  the  towers  are  separated  from  the  cable  and  each 
piece  is  written  out  as  a separate  feature.  A wide  river  is  represented 
by  two  separate  features,  one  defining  its  right  bank  and  another  defining 
its  left  bank.  In  addition  to  maintaining  the  two  banks,  RSS1  combines 
them  into  a third  feature  which  is  a closed  polygon  representing  the  water 
sur face. 


1 • 3 Output 

RSS1  outputs  the  processed  planimetry  file  to  disc.  Each 
feature  is  represented  by  a record  with  the  following  format: 

Word  1:  Feature  code 

Word  2:  Total  number  of  coordinates  in  the  record  i.e.  twice 

the  number  of  points  N. 

Word  3-word  2N-I-2:  X,Y  data 

RSS1  also  outputs  a printout  identifying  the  type  and  location  of 

all  features  on  the  map  sheet.  A sample  of  this  printout  is  given  in 
Figure  II. 2. 


1.4  Timing  and  Cost 

For  a typical  map  sheet  containing  131  features,  RSS1 
requires  92  cp  seconds  execution  time  and  3*s  minutes  wall-clock  time 
on  ETL's  CDC  6400.  Run  cost  is  approximately  twenty-four  dollars. 

2.  Program  RSS2 

2.1  Input 

Program  F.SS2  requires  two  inputs.  The  first  is  the  output 
file  from  RSS1.  This  transfer  is  accomplished  through  the  use  of  SCOPE 
control  cards  and  will  be  discussed  in  Chapter  III.  The  second  is  a single 
card  input.  The  format  for  this  card  is  120,  15X,  15,  8X,  12  and  its  con- 
tent is  as  follows: 
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Columns  1-20: 


contain  the  map  scale,  e.g.  if  the  map  scale  is 
1:100,000,  this  field  contains  the  number  100,000. 

Columns  21-35:  blank 


Columns  36-40:  contain  the  size  of  the  grid  resolution  elements  used 

to  label  points  in  a region  of  planimetry  data.  This 
value  must  be  expressed  in  units  of  meters  x 1000. 

Columns  41-48:  blank 


Columns  49-50:  contain  the  region  size,  i.e.  the  number  of  resolution 

elements  along  a region  edge.  The  restrictions  on  this 
value  are  (1)  it  must  not  exceed  48  and  (2)  it  must  be 
a multiple  of  4. 


2.2  Description  of  Processing 

As  mentioned  previously,  RSS2  generates  the  data  strips  used 
to  describe  planimetry  features  in  RSS.  The  process  of  strip  genera- 
tion can  be  roughly  described  as  follows. 

First,  the  map  scale  and  grid  size  are  used  to  express  all  X,Y  data 
pairs  in  units  of  the  grid  size.  Now  consider  two  consecutive  data  points 
(Xl'Yl>  and  (X^.Y^).  To  construct  the  data  strips  corresponding  to  these 

two  points  begin  by  calculating  the  inverse  slope. 


S = 


Ax 

zj  Y 


VWY 


2 


The  number  of  strips  to  be  generated  is  given  by  Y1-Y2  • i-e.  each  strip 

is  of  unit  width.  The  first  strip  begins  at  Y=Y1»  X=X^  and  is  of 

length  X+S,  since  X changes  by  a distance  equal  to  S when  Y changes  by  1. 
The  second  strip  has  Y=Y,-1,  X=X1+S+1  and  ^X=S,  and  so  on  for  the  |Y1_Y2| 

strips  generated  between  the  two  points.  The  procedure  is  then  repeated 
for  (X2 .Y2)  and  CX3, Y3> , etc. 

This  description  is  of  course  an  oversimplification.  The  process 
actually  does  not  employ  integer  arithmetic  since  S may  not  be  an  inte- 
ger. Also,  the  sign  of  the  slope  is  a pertinent  quantity  since  it  deter- 
mines whether  successive  strips  move  to  the  left  or  to  the  right.  However, 
the  essence  of  the  method  is  as  described  above. 

The  data  strips  for  each  feature  are  sorted  by  Y and  then  x.  The 
placement  of  the  data  in  this  form  is  required  by  RSS5  which  merges  the 
strips  for  overlapping  features. 


2.3 


Output 


RSS2  outputs  two  disc  files.  The  first  file  outputs  the 
data  for  each  feature  in  blocks  of  1800  words.  The  first  six  words  of 
the  first  block  contain  general  information  about  the  feature,  with  word 
2 containing  the  feature  code.  Beginning  in  word  7 are  the  values  of  the 
quantities  describing  the  data  strips.  The  format  of  this  information  is 
as  follows: 


WORD  CONTENT 

6+4N-3  Y-value  for  the  Nth  strip 


6+4N-2 

6+4N-1 


6+4N 


X-value  for  the  beginning  of  the  Nth  strip 

Ax  value  for  the  Nth  strip;  equals  the 
total  number  of  resolution  elements  covered 
by  the  strip.  The  endpoint  of  the  strip  is 
then  at  X+/^X-1 

Specularity  angle  for  the  Nth  strip.  This 
quantity  is  used  to  determine  abnormal  radar 
return  qualities  and  is  not  presently  used 
in  RSS,  although  it  is  calculated. 


The  second  file  is  a two-word  file  containing  the  following  map  size 
values : 


Word  1 Grid  Resolution  element  size  (meters  X 1000) 

Word  2 Region  size 

RSS2  also  outputs  a printout  of  the  header  information  for  each  feature. 

2.4  Timing  and  Cost 

Because  of  the  large  processing  requirements,  RSS2  is  an  expen- 
sive and  time-consuming  program  to  operate.  For  the  file  of  190  features  1465  cp 
seconds  and  42  minutes  wall  clock  time  were  required  for  execution.  Run  cost  was 
353  dollars. 

3.  Program  RSS 3 

3.1  Inputs 

Program  RSS3  requires  as  input  the  output  files  from  RSS2. 

No  card  input  is  needed.  Information  linkage  is  accomplished  through  the 
use  of  SCOPE  control  cards  as  described  in  Chapter  III. 

3 # 2 Description  of  Processing 

RSS3  assigns  the  planimetry  data  strips  to  regions.  Each 
planimetry  region  is  a square  whose  sides  represent  a ground  distance  equal 
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to  the  product  of  the  region  size  and  the  resolution  grid  size.  Numbering 
of  these  regions  begins  in  the  southwest  corner  of  the  map  with  the  assign- 
ments being  made  sequentially  left  to  right  and  bottom  to  top. 


The  procedure  used  to  assign  the  planimetry  strips  to  their  respective 
regions  is  a simple  one.  Consider  for  example,  a strip  with  the  following 
characteristics : 


Region  size  = 32 
Resolution  = 156.25  feet 


This  means  that  the  strip  begins  117,187  feet  north  (750  X 156.25) 
and  136,718  feet  east  (875  x 156.25)  of  the  map  origin.  We  have  assumed 
that  each  region  contains  32  increments  along  each  edge  so: 


750/3 


23.43750 


If  we  assume  that  each  horizontal  row  of  regions  can  contain  at  most  331 
regions,  then  the  region  number  for  the  beginning  of  the  strip  is: 


(23)  + 28  = 7641 


The  (X,Y)  coordinates  of  the  beginning  point  relative  to  the  first  qrid 
element  in  the  region  is: 


Y = 0.43750  * (32)  = 14 
X = 0.34375  * (32)  = 11 


Since  ;\X  = 24,  the  end  point  of  the  strip  is  at  X = 11+24-1  = 34.  The 
maximum  address  within  a region  is  32,  so  this  strip  overflows  into  the 
next  region  to  the  east.  We  therefore  wind  up  with  two  strips  as  follows 


7641 

14 

11 

22 


7642 

14 


RSS3  also  assigns  merge  priorities  for  each  feature  code.  This  is 
done  through  the  use  of  a translation  table  which  is  hard-coded  as  a DATA 
statement. 


The  output  from  RSS3  is  placed  on  disc  in  the  form  of  card 
Each  image  is  labeled  by  the  region  number  and  feature  code.  It 


contains  the  merge  priority  for  the  feature  and  the  descriptive  information 
for  up  to  seven  planimetry  strips. 


The  only  printed  output  from  RSS3  ls  a message  indicating  the  end  of 
processing. 

3.4  SORT  Routine  (RSS4) 


Program  RSS3  employs  the  CDC  SORT/MERGE  package  to  order 
the  card  images.  For  the  purposes  of  this  document,  we  will  call  this 
routine  RSS4. 


In  preparation  for  RSS5,  the  card  images  are  sorted  in  the  following 


order: 


ORDER 


ATTRIBUTE 


1 

2 

3 

4 


Region  Number 
Merge  Priority 
Y-value  of  first  strip 
X-value  of  first  strip 


Although  the  reader  is  directed  to  the  CDC  SORT/MERGE  manual  for  com- 
plete details  on  the  operation  of  the  SORT/MEP.GE  package,  we  present  here  a 
listing  and  brief  description  of  RSS4. 

SORT,  VAR=DISC 
BYTES I ZE , 6 

FILE, SORT=TAPE 1 , OUTPUT=TAPE  2 

FIELD, REGION (1.1, 5, DISPLAY) . PRIORITY (6 . 1 , 2 , DISPLAY) , 

, INTY (13.1,2, DISPLAY)  , INTX (15.1,2, DISPLAY ) 

KEY, REGION (A, OWN ) , PRIORITY < A, OWN)  ,INTY(A,OWN)  ,INTX(Z,OWN) 
SEQUENCE, OWN (#,0,1, 2, 3, 4, 5, 6, 7, 8, 9) 

OPTIONS, RETAIN 
END 


• a 


The  FIELD  card  lists  the  attributes  to  be  sorted  in  order  of  decreasing 
priority.  It  indicates  the  format  of  the  data,  which  in  our  case  is  CDC 
DISPLAY  code,  and  its  location  within  the  record.  The  KEY  card  gives  the 
direction  of  the  sort  (i.e.  the  "A"  indicates  that  the  numbers  are  to  be 
placed  in  ascending  order)  and  the  ordering  scheme,  which  in  our  case  is  the 
sequence  OWN.  The  SEQUENCE  card  is  used  to  define  an  ordering  sequence 
other  than  one  of  the  standard  sequences.  This  card  is  required  in  RSS4 
because  a blank  in  display  code  is  represented  by  55B  and  this  value  is 
greater  than  that  assigned  to  the  integers.  Without  a special  sequence, 
a number  like  b#1023  would,  for  example,  be  placed  after  kS16486  in  an 
ascending  order  sort  of  the  card  images. 

Clearly,  the  output  from  RSS4  is  in  the  same  format  as  that  from 
RSS3.  Only  the  ordering  of  the  records  is  changed. 
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3.5 


Timing  and  Cost 


(RSS3  and  RSS4) 


For  our  typical  data  base  of  190  features  the  running  time 
for  the  combination  program  is  231  cp  seconds  and  2 minutes  wall  clock  time. 
Cost  for  the  run  is  thirty-four  dollars. 


4 . Program  RSS5 


4. 1 Inputs 


RSS5  requires  the  output  file  from  RI.MS4  and  the  two  word 
disc  file  from  RSS2.  It  does  not  require  any  card  input. 

4. 2 Description  of  Processing 

This  program  merges  the  planimetry  strips  so  that  no  strips 
overlap.  The  procedure  has  been  discussed  in  Chapter  I and  is  conceptually 
quite  simple.  The  records  containing  the  planimetry  strips  have  already 
been  sorted  by  region  number  and  merge  priority.  The  merge  priority  scheme 
is  set  up  so  that  those  features  which  may  be  overwritten  appear  before 
those  which  may  not.  For  example,  for  a given  region,  strips  describing 
a lake  appear  before  those  describing  an  island,  and  those  describing  a 
road  occur  after  those  describing  the  island.  The  scheme  is  presently 
imperfect  (consider  an  island  with  a lake  on  it)  but  can  be  expected  to 
work  in  the  vast  majority  of  cases. 

In  essence,  RSS5  simply  copies  the  planimetry  strips  into  a core 
array  in  the  order  in  which  they  appear  on  the  input  file.  This  array 
becomes  an  image  of  the  region  structure  as  successive  strips  are  copied 
into  it.  An  empty  record  is  written  for  those  regions  which  contain  no 
planimetry  data. 

As  an  example  of  the  operation  of  RSS5,  consider  two  planimetry 
strips,  one  being  part  of  a lake  and  the  other  being  part  of  a city. 

Take  the  location  and  size  of  these  strips  to  be  as  follows: 


Lake  Strip 

Y = 25 
X = 11 
LlX  = 6 


City  Strip 

Y = 25 
X = 8 
X = 15 


Clearly,  this  portion  of  the  lake  overlaps  the  city.  RSS5  takes 
these  strips  and  generates  three  (3)  strips,  two  belonging  to  the  city 
and  one  belonging  to  the  lake.  The  location  and  size  of  these  strips 
is  as  follows: 


Lake  Strip 


City  Strip 


Y 

X 

X 


25 

11 

6 


Notice  that  these  strips  do  not  overlap. 


Y = 25 
X = 17 
X = 6 
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4.3 


Output 

As  just  mentioned,  RSSS  outputs  a record  for  each  region  on 
the  map  regardless  of  whether  or  not  that  region  actually  contains  any  plani- 
metry data.  The  first  word  of  the  record  contains  the  total  number  of  strips 
contained  in  the  region,  while  the  remaining  words  contain  the  information 
describing  the  strips.  Each  strip  is  described  by  one  word,  with  the  data 
being  packed  as  follows: 

BITS  INFORMATION 

31-36  Y-coordinate  of  strip 

25-30  X-coordinate  of  beginning  of  strip 

19-24  . X for  strip 

7-18  Specularity  angle  for  strip 

1-6  Feature  code  for  strip 

The  blocks  are  written  to  a random  access  file  with  the  location  of  each 
block  being  labeled  by  the  corresponding  region  number. 

The  printed  output  from  RSS5  includes  a printout  of  all  regions 
that  were  padded  with  empty  data  to  yield  a region  size  with  a perfect 
square,  and  an  end-of-processing  message. 

Timing  and  Cost 

This  program  required  216  cp  seconds  and  2 minutes  wall  deck 
time  to  process  the  strips  for  190  features.  Run  cost  was  forty-four  dollars. 

5.  Program  RSS6 

5.1  Inputs 

The  input  to  this  program  is  a tape  containing  elevation  values 
for  the  area  covered  by  the  map  and  the  two  word  disc  file  output  from  RSS2. 

The  first  record  on  the  tape  is  a 36  word  header.  Each  entry  in  the 
header  is  an  integer  representation  of  the  actual  data  obtained  by  multi- 
plying the  value  by  1000  and  rounding.  The  first  20  words  contain  infor- 
mation pertaining  to  the  location  on  the  earth  of  the  area  from  which  the 
data  was  obtained  and  the  values  of  parameters  required  in  producing  the 
cartesian  projection  of  the  map.  These  quantities  will  not  concern  us  here. 
The  last  16  entries  are  as  follows: 

Word  21:  The  total  number  of  records  (profiles)  on  the  tape, 

not  counting  the  header. 

Word  22:  The  total  number  of  elevation  points  per  record. 

Elevations  are  unsigned  15  bit  integers,  in  feet, 
packed  4 per  word. 

Word  23-25: 


Geographic  latitude  in  degrees,  minutes,  and  seconds 
of  the  radar  target. 


m 


u 

0 

D 

D 

D 

Q 

13 

D 

n 

D 

D 

D 

D 

D 

D 

n 

n 

n 


Word  26:  Ground  distance  in  meters  between  the  UNAMACE  point 

of  tangency  north  to  tne  target. 

Word  27-29: 

The  longitude  in  degrees,  minutes,  and  seconds  of  the 
radar  target. 

Word  30:  The  ground  distance  in  meters  from  the  UNAMACE  point 

of  tangency  west  to  the  target. 

Word  31:  The  ground  distance  in  meters  from  the  first  elevation 

profile  (upper  NW  corner  of  the  elevation  data)  north 
to  the  radar  target.  A negative  value  indicates  that 
the  target  is  south  of  the  first  elevation  point. 

Word  32:  The  ground  distance  in  meters  from  the  first  elevation 

profile  west  to  the  target.  This  number  is  also  negative. 

Word  33:  Target  ID  Code 

Word  34:  Ground  spacing  between  points  of  a profile  in  meters. 

Presently  this  value  must  be  equal  to  the  grid  resolu- 
tion size  output  by  RSS2. 

Word  35:  Ground  spacing  between  profiles  in  meters.  Equals 

word  34  at  present  time. 

Word  36:  Zero  fill  word. 

Following  this  header  are  a series  of  records  - one  for  each  profile  line. 
Each  profile  contains  the  elevation  values  for  a west-to-east-running 
strip  of  terrain,  with  the  first  profile  corresponding  to  the  northern- 
most part  of  the  map  (a  contrast  to  the  SW  origin  for  digitizing  the  plani- 
metry) . 


Description  of  Processinc 


The  input  data  base  is  characterized  by  resolution  elements 
equal  to  those  of  the  planimetry  file  with  the  data  being  arranged  in  hori- 
zontal strips  covering  the  entire  width  (east-west  dimension)  of  the  map 
sheet,  RSS6  reformats  the  data  in  these  strips  into  regions  with  dimen- 
sions equal  to  the  dimensions  used  by  the  planimetry  file.  The  regions  are 
numbered  from  the  lower  left  corner  of  the  map  just  as  in  the  planimetry 
data  file. 


For  example,  if  the  region  size  is  defined  as  24  by  RSS2,  then  24 
profiles  are  read  at  a time  into  a core  array.  From  word  22  of  the  header 
it  can  be  determined  how  many  regions  (in  the  horizontal  direction)  are 
defined  by  this  data.  Similarly,  word  21  of  the  header  permits  the  deter- 
mination of  the  Y-coordinate  (i.e.  number  of  regions  in  the  Y-direction) 
of  each  row.  From  this  information  it  is  possible  to  assign  a region  number 
to  each  group  of  24  X 24  data  points  thereby  defining  a grid. 
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5.3  Output 

The  information  for  each  region  is  stored  in  a one-dimensional 
array.  The  first  word  contains  the  region  number  while  the  remaining  words 
contain  the  elevation  data  packed  four  per  word.  The  packing  order  is  such 
that  if  there  are  N resolution  elements  along  a region  edge,  each  sequential 
group  of  N/4  words  contains  the  elevation  data  for  one  row  of  resolution  ele- 
ments within  the  region  beginning  at  the  western  edge,  with  the  southern-most 
row  appearing  first,  i.e.: 

WORD  1 REGION  NUMBER 

Word  N+l  Elevation  value  for  row=(N-l)/IK  +1 

column  = N-IR  * (row-1) , where  IR  is 
the  number  of  resolution  elements 
along  a region  edge. 

During  processing,  each  region  record  is  placed  on  a random  access  file.  When 
all  regions  have  been  processed,  the  records  are  read  back  into  core  in  region 
order  and  rewritten  to  the  disc. 

5.4  Timing  and  Cost 


Program  RSS6  requires  117  cp  seconds  and  16  minutes  wall 
clock  time  for  execution.  Run  cost  is  approximately  forty-three  dollars. 

6,  Program  RSS7 


6.1  Input 

Program  RSS7  requires  tape,  disc,  and  card  input. 

The  tape  required  is  the  output  tape  from  RSS6,  while  the 
disc  file  is  the  output  of  RSS5  and  the  two  word  file  from  RSS2.  The 
card  input  supplies  certain  NAMELIST  quantities  and  their  meanings  are  as 
follows : 


QUANTITY 


MEANING 


IFREQ 


The  number  of  radial  scan  lines  to  be  constructed 
per  degree.  Fcr  example  IFREQ=2  will  result  in 
720  scan  lines  separated  by  1/2°. 

Default:  IFREQ=2 


ITERR 


Control  variable.  INTERR=1  causes  the  planimetry 
data  to  be  merged  with  the  terrain  data.  INTERR=0 
causes  the  terrain  data  to  be  ignored. 

Default:  INTERR=1 
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ALTTDE 


6.2 


The  maximum  altitude  of  the  radar  above 

the  radar  target.  The  maximum  value  of  this  field 

is  32000 

Default:  ALTTDE=  32000  feet 

X-coordinate  of  the  radar  target  in  meters  relative  to 
the  southwest  corner  of  the  map 

Y-coordinate  of  the  radar  target  in  meters  relative  to 
the  southwest  corner  of  the  map 

Description  of  Processing 


The  purpose  of  RSS7  is  to  merge  the  planimetry  data  file 
from  RSS5  with  the  terrain  data  file  from  RSS6  and  to  determine  certain 
map  related  quantities  needed  to  construct  the  radar  scene. 

Using  information  from  the  terrain  file  header  record,  the  region 
address  of  the  radar  target  is  determined.  With  this  point  as  its  center, 
an  imaginary  circle  is  then  constructed  whose  radius  is  equal  to  the  ground 
range  of  the  radar.  A square  circumscribed  about  this  circle  will  then  con- 
tain all  of  the  regions  required  to  produce  the  radar  scene  at  the  chosen 
altitude;  the  rest  may  be  discarded.  This  procedure  is  illustrated  in 
Figure  II. 3.  For  each  region  overlapped  by  the  circumscribed  square,  the 
terrain  and  planimetry  data  are  merged  into  a single  record  labeled  by  the 
region  number. 
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Output 


RSS7  outputs  two  files.  The  first  contains  two  five-word 
records  containing  descriptive  information  required  in  the  reconstruction 
of  the  radar  scene. 


RE CORD 1 

WORDl 

W0RD2 

WORD  3 

WORD4 

WORD5 

RECORD 2 
WORDl 


CONTENTS 

Ground  range  of  the  radar  in  meters. 

Radar  height  (in  meters)  above  target  ground 
altitude 

X-coordinate  of  the  radar  target  in  meters  rela- 
tive to  the  southwest  corner  of  the  nap. 

Y-coordinate  of  the  radar  target  in  meters  rela- 
tive to  the  southwest  corner  of  the  map. 

Elevation  (in  meters)  of  the  radar  target  above 
sea  level.  The  sum  of  this  value  and  WORD2 
gives  the  radar  altitude  above  sea-level. 


Maximum  number  of  regions  in  the  X-direction  of 
the  map  sheet;  presently  set  at  331. 
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A 


WORD2 

The 

WORD  3 

The 

WORD4 

The 

Y 

WORDS 

The 

resolution  unit  size  in  meters  X 1000. 
value  of  I FREQ 

region  size  dimension  in  both  the  X and 
direction. 

value  of  ITERR 


The  second  file  is  a random-access  file  containing  the  merged  plani- 
metry-terrain data  records.  These  records  are  labeled  by  region  number  and 
have  the  following  format. 


Region  number 

Descriptive  information  for  the  N strips  contained 
in  the  region.  Each  word  contains  the  information 
for  one  strip,  as  outputted  by  RSS5. 

Packed  elevation  data  (four  per  word)  as  outputted 
from  RSS7.  The  ordering  of  this  elevation  data, 
and  its  relation  to  the  corresponding  grid  loca- 
tions within  the  region  is  the  same  as  that  of 
the  output  from  RSS6. 

Figure  II. 4 shows  a sample  printout  from  RS57. 

4 Timing  and  Cost 

Program  RSS7  executes  in  approximately  36  cp  seconds  with 
IFREQ=2.  Run  cost  is  eleven  dollars. 

7 . Program  RSS8 

7 • 1 Input 

Program  RSS8  requires  the  two  output  files  from  rss7 
the  two  word  disc  file  from  RSS2. 


W0RD1 

W0RD2 

WORD  N+l  - 
WORD  N+IIC 


7 . 2 Description  of  Processing 

This  program  creates  the  radar  scan  lines  from  which  the 
final  radar  scene  is  constructed.  The  procedure  for  doing  this  will  now 
be  outlined. 


In  polar  coordinates,  any  point  on  the  map  may  be  described  by  the 
following  expression: 

X = X + r cose 

p t p 

Y = Y + r SIN6 
p t p 

where  (X  ,Y  ) is  the  location  of  the  target,  r the  radial  distance  from 
t t p 
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the  target,  and  0 the  polar  angle.  The  quantity  IFREQ  defines  the  values 
of  6 to  be  used  in  the  generation  of  the  scan  lines.  These  values  are: 

0^  = n/IFREQ , 0 ^n^360*IFREQ 

and  result  in  the  generation  of  a circle  consisting  of  360*IFREQ  scan 
lines  separated  by  an  angular  increment  of  1/IFREQ  degrees.  The  radar 
range  R is  expressed  in  units  equal  to  the  grid  spacing  for  the  plani- 
metry data.  The  variable  r therefore  has  the  range  0 <r  ^R. 

P P 

With  the  parameters  thusly  defined,  the  generation  of  the  sweep 
lines  proceeds  as  follows.  We  asusme  that  the  initial  value  of  0 =0, 

i.e.,  n=0  and  that  we  begin  each  scan  line  r = 1. 

P 

A)  Calculate  COS0  and  SIN8 

n n 

B)  Calculate  X and  Y 

P P 

C)  Calculate  the  region  in  which  (X  ,Y  ) lies 

P P 

D)  If  the  record  for  this  region  is  already  in  core, 
go  to  (G) 

E)  Read  the  record  for  this  region  into  core 

F)  Copy  the  planimetry  strip  information  for  this  region 
into  a IK  X IK  core  array  which  serves  as  an  image  of 
the  region.  Radar  return  intensity  levels  are  assigned 
in  correspondence  with  the  feature  codes  by  using  an 
internal  translation  table.  The  radar  return  intensity 
for  resolution  elements  not  covered  by  planimetry  data 
is  set  equal  to  a predetermined  background  intensity. 

The  result  is  that  ARRAY  (l,j)  contains  the  radar  re- 
turn intensity  (from  planimetry  only)  for  the  resolution 
element  at  X=I , Y=J. 

G)  Calculate  which  resolution  element  contains  (X  ,Y  ) . 

P P 

Denote  this  by  ARRAY  (I  ,J  ) 

P P 

H)  Copy  ARRAY  (I  ,J  ) and  the  corresponding  terrain  value 

P P 

from  the  input  region  record  into  a linear  array  which 
stores  the  information  for  the  scan  line. 


I)  Let  r 


J)  If  r ^R  go  to  (B. ) 

K)  Let  0 — >0  + 1/IFREQ 

L)  If  0>  360°  we  are  finished.  If  & ^360°  go  to  (A.) 
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FIGURE  II  .4  - SAMPLE  F r UT^’T  PROM  RSS7_ 
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FIGURE  1 1. 5 - SAMPLE  PRINTOUT  FROM  RSS8 
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Output 


The  output  file  from  RSS8  consists  of  360*1FREQ  records  of 
4000  words  each.  The  file  is  sequential  in  nature  with  the  records  being 
ordered  by  their  corresponding  angle  variable  9 i.e.  the  record  for  the 
scan  line  at  9 = 0 is  the  first  one  on  the  file  while  that  for  9=360  is 
the  last  one. 

The  first  2000  words  contain  the  intensity-of-return  values  for  each 
4,  0 <r^R.  Clearly,  if  R <2000  the  remaining  words  are  blank.  Words  num- 
bered 2000  + I,  1^  2000,  contain  the  corresponding  terrain  elevation 

values.  For  example,  if  IFREQ=2,  then  the  fourth  record  has  the  following 
significance. 

A)  It  corresponds  to  9 = 4/IFREQ=2° 

B)  The  Ith  word  and  the  I+2000th  word  respectively  contain 
the  strength  of  return  and  the  terrain  elevation  for 
the  point: 

X = X + I *C0S2° 
p T 

y = Y + I*SIN2° 
p T 

The  printout  from  RSS8  is  shown  in  Figure  II. 5 and  is  self-explanatory. 
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Timing  and  Cost 


Program  RSS8  executes  in  approximately  1796  cp  seconds.  Run 
cost  is  approximately  390  dollars. 

8-  Program  RSS9 

8.1  Input 

Program  RSS9  requires  the  output  file  from  RSS8  and  the 
parameter  file  from  RSS7,  the  two  word  disc  file  from  RSS2,  and  a one 
card  input  containing  the  radar  altitude. 


8.2 


Description  of  Processing 


RSS9  processes  the  data  base  to  take  into  account  certain 
geometrical  effects  required  to  make  the  final  image  similar  to  that  which 
may  be  expected  from  a real  radar.  Prior  to  RSS9  processing,  the  only 
radar  return  strengths  included  in  the  data  base  are  certain  intrinsic  re- 
turns assigned  to  the  planimetry  strips  on  the  basis  of  feature  code  in 
RSS8.  AH  of  the  background  ("between"  the  scan  lines)  has  been  assigned 
the  same  return  intensity  and  this  is  clearly  inaccurate  because  it  fails 
to  take  into  account  the  ground  topography.  RSS9  assigns  a radar  return 
strength  to  every  resolution  element  within  20  nautical  miles  of  the  target 
location  based  on  the  detailed  structure  of  the  earth's  surface  in  the  area 
under  consideration. 
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Three  specific  effects  are  considered.  The  first  is  the  Lambert's 
Law  correction  and  is  illustrated  in  Figure  II. 6.  The  mathematical  state- 
ment of  this  law  is  that  the  return  from  any  given  point  on  the  ground  is 

reduced  by  a factor  COS0  where  0 is  the  angle  between  the  incident  radar 

L L 

beam  and  the  normal  to  the  terrain  at  the  point  of  interest.  This  is  to 
say,  if  I is  the  intrinsic  intensity  of  the  background  or  of  the  planimetry 
feature  located  at  a given  point,  then  the  actual  radar  return  intensity  is 
given  by: 


I = 
r 


COS0. 


The  implementation  of  this  effect  in  RSS9  is  quite  straightforward. 

For  each  value  of  r,  the  slope  of  the  line  through  the  points  r-1  and  r+1 
is  calculated  from  the  corresponding  terrain  elevation  values  thus  yielding 
the  angle  0^.  The  vector  perpendicular  to  this  line  defines  the  normal  to 

the  surface  at  r.  The  declination  angle  of  the  radar  beam  0 can  be  ealcu- 

D 

lated  from  r,  the  terrain  elevation  at  r and  the  radar  altitude.  Using 

these  angles,  the  angle  0 can  be  determined  from  elementary  trigonometry. 

L 

A second  effect  taken  into  account  is  that  of  shadowing.  The  pertinent 
geometric  considerations  are  illustrated  in  Figure  II. 7.  The  nature  of  the 
effect  is  simply  that  points  on  the  ground  may  be  in  the  shadow  of  terrain 
peaks  lying  between  the  radar  and  the  ground  point  of  interest.  Such  points 
will  be  invisible  to  the  radar  and  must  therefore  be  assigned  a zero  intensity 
of- return. 


The  third  radar  effect  that  RSS9  handles  is  the  altitude  layover  effect 
This  effect  shifts  the  location  of  the  terrain.  The  length  and  direction  of 
the  shifting  of  each  processed  point  of  the  terrain  depends  on  the  height  of 
that  point.  This  effect  is  implemented  by  considering  each  point  R along  a 
radial  line  and  computing  a new  altitude  layover  point  Rl.  Given  (1)  the 
distance  of  the  point  (R)  away  from  the  target  point  along  a radial  line, 

(2)  the  height  of  the  radar  above  sea  level  (A) , and  (3)  the  elevation  of 
the  point  R above  sea  level  (E) ; the  distance  (D)  from  the  radar  to  the 
elevated  point  along  the  radial  can  be  computed  as  follows: 


n (A- 


E)2  + 


D 


Using  a principle  of  elementary  trigonometry,  we  know  that  a tangent 
line  drawn  the  length  of  the  elevation  of  point  R to  some  point  Ri  at  target 
level  alo-,g  the  radial  line  yields  a distance  D1  equal  in  length  to  distance 
D.  The  point  Rl  can  be  computed  similarly  to  the  above  equation.  Having 
computed  D=Dl , and  given  the  radar  altitude  above  target  level  (ALT),  the 
point  Rl  computed  at  target  level  is  as  follows: 

f 2 2 2 

= v ALT  + (A-E)  + R 


Rl 
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FIGURE  1 1. 6 - DEFINITIONS  OF  LAMBERT'S  LAW  QUANITIES 


FIGURE 


RADAR  GROUND  COVERAGE  AS  A FUNCTION  OF  ALTITUDE  AND  APERTURE 


In  addition  to  the  three  radar  effects  listed  above,  RSS9  also  scales 
the  scan  lines  to  insure  that  all  generated  radar  scenes  are  the  same  size. 
The  need  for  such  an  operation  is  illustrated  in  Figure  II. 8.,  and  is  due  to 
the  fact  that  with  a fixed  radar  aperture  the  ground  distance  covered  by  the 
radar  will  be  the  function  of  the  altitude  of  the  radar.  Since  we  will  asso- 
ciate one  resolution  element  with  one  pixel  on  the  DICOMED  display,  a smaller 
ground  range  will  result  in  a progressively  smaller  display.  Subroutine  SIZE 
within  RSS9  therefore  expands  or  condenses  the  image  by  a factor  equal  t.o: 

F = 1000/NQ.  OF  RESOLUTION  ELEMENTS  IN  THE  GROUND  RANGE 

where  we  have  chosen  1000  resolution  elements  (which  is  equivalent  to  19.194 
nautical  miles  for  a radar  located  32000  feet  above  ground  level)  as  our 
desired  image  size.  This  is  done  prior  to  processing  by  subroutine  LAMBERT 
which  includes  the  radar  effects  listed  above.  SIZE  simply  assigns  the  ele- 
vation and  radar  return  for  distance  r to  distance  F*r.  For  points  between 

F*r  and  F*(r+1),  return  intensities  are  set  equal  to  the  background  value 

and  elevations  are  computed  by  linear  interpolation. 

8 . 3 Output 

The  radial  sweep-line  format  of  the  data  in  RSS9  is  quite 
suited  to  the  task  of  incorporating  radar  effects  since  such  effects  are 
functions  of  the  radial  distance.  However,  for  purposes  of  display  on  the 

DICOMED  plotter,  the  image  information  must  be  in  a raster  format.  As  a 

preliminary  to  this  raster  conversion,  RSS9  outputs  the  data  as  a large 
number  of  small  records,  each  corresponding  to  a single  point  on  the  image. 
This  is  done  in  subroutine  REFMT.  Each  record  contains  the  following  infor- 
mation : 


W0RD1 


Number  containing  encoded  (X,Y)  coordinate 
information 


WORD2 


Intensity  value  from  0-63.  Here  intensity 
zero  corresponds  to  maximum  radar  return  (white) 
while  a 63  corresponds  to  zero  return.  This 
system  is  the  same  as  used  on  the  DICOMED  plotter. 


This  data  is  written  to  four  files,  each  of  which  contains  approximately 
1/4  of  the  image.  This  is  done  to  speed  up  the  sorting  of  the  points  required 
to  order  them  for  raster  plotting. 

The  printout  from  RSS9  is  shown  in  Figure  II. 10.  The  input  statistics 
indicate  the  distribution  of  color  codes  as  outputted  by  RSS8.  In  the  case 
presented  here,  a color  code  of  2 was  used  for  the  background. 

8.4  Timing  and  Cost 


RSS9  requires  approximately  1009  cp  seconds  and  6 minutes 
wall  clock  time  for  execution.  Run  cost  is  approximately  133  dollars. 


rfLMSS  SUCCESSFUL  lN0,  OISPLAV  rh[ 

INPUT  statistics 

133U75  PIXTLS  WIIH  uOLOR  COOT.  ;• 

649  PIXFLS  WITH  COLOR  COOr  j t 
lo<C  PULLS  WITH  COLOR  Cujt  l,; 

JO'.  RUTLS  WIIH  CJluk  C III.  l, 

3UTUUT  STATISTIC., 

p PIXILS  WITH  COLOR  CM'!  ft 
lb  PULLS  WITH  COLOR  C'jJF  0 
2 9 PIXILS  HUH  COLO1-'  CUUF  10 
26  PIXLLS  WITH  COLOR  C)Jt  11 
60  PIXELS  WITH  LOLUE  CUJL  1? 

32  PIXLLS  WITH  CULuP  COIjr  13 
40  PIXLLS  WITH  COLOR  CUUL  14 
6 2 PIXLLS  WITH  COLOR  COJE  15 
fi  PIXLLS  WITH  COLOP  COOL  lb 

80  pixfls  with  color  cool  it 

85  PIXlLS  WITH  COLOR  COOE  18 
9b  PIXELS  WITH  COLOR  COOE  19 
94  PIXELS  WITH  COLOR  COOE  20 
149  P I Xc.L  S WITH  COLOR  COOr  21 
lob  PIXtLS  WITH  COLUR  COOE  22 
18b  PIXELS  WITH  COLOR  COOE  23 
239  PIXLLS  WITH  COLUR  CJUE  24 
279  PIXELS  WITH  oOLOR  CODE  25 
Jib  PIXELS  WITH  CULUk  COOE  2b 
389  PIXELS  WITH  COLOR  COOE  27 
51)2  3 IXELS  WITH  COLOR  COJF  28 
558  PIXELS  WITH  COLOR  COOE  29 
087  PIXELS  WITH  COLOR  CJOE  30 
9 7 J PIXELS  WITH  COLOR  COUE  31 
1407  PIXELS  WITH  COLOR  COOE  J? 

232 b PIXELS  WITH  COLOR  COOF  33 
369b  PIXlLS  WITH  COLOR  COOF  34 
4849  PIXLLS  WITH  COLOR  COOE  35 
5932  PIXFLS  WITH  COLOR  COOE  36 
9322  p I XEL h WITH  COLOR  COOE  37 
13391  PIXELS  WITH  COLOR  COOE  38 
21559  PIXtLS  WITH  CCLOk  CODE  39 
30185  PIXELS  WITH  COLOR  COJF  40 
15275  PIXELS  WITH  COLOR  COJE  41 
7029  PIXELS  WITH  COLOR  COJE  42 
3936  PIXELS  WITH  oOLJP  C0Oc  43 
2418  PIXELS  WITH  COLOR  COJE  44 
1 l 32  PIXELS  WITH  COLOR  COOE  45 
1029  PIXELS  WITH  COLOR  COO;  4f 
790  PIXELS  WITH  COLO'-’  C9Jr  47 
594  PIXELS  WITH  COLOR  COJE  48 
437  P I XCL5  WITH  COLOR  OJJE  49 
359  ’IXELS  WITH  CULUR  COJE  50 
331  PIXELS  WIIH  COLOR  COOE  51 
308  PIXELS  WITH  COLOR  CJOE  52 
569  PIXELS  WITH  COLOR  COOE  53 
383  P IXELS  WITH  COLOR  COJF  54 
207  PIXlLS  WITH  COLOP  COJE  55 
235  PIXLLS  WITH  COLOR  CJOE  56 
161  PIXELS  WIIH  COLOR  CJOE  57 
14>  PIXELS  WIIH  COLOR  COJF  58 
102  PIXtLS  WITH  COLOR  COOE  59 

105  PIXFLS  WITH  COLOR  COJc  61 

106  PIXELS  WITH  COLOR  COJF  61 
441  PIXELS  WITH  COLOR  C9JF  o2 
970  PIXELS  WIIH  COLOR  COJF  63 

924  PIXELS  WERE  ASSTCNtJ  TO  T RE  SHAOOW 


FIGURE  II. 10-  SAMPLE  PRINTOUT  FROM  RSS9 
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9.  Program  SORT 


This  program  uses  the  CDC  SORT/MERGE  package  to  combine  the  four 
output  files  from  RSS9  into  a single  file  ordered  as  follows:  the  records 

are  sorted  first  by  Y and  then  by  X.  Because  the  X,Y  information  for  each 
point  is  encoded  into  a single  number,  this  can  be  accomplished  via  a single 
level  sort.  The  procedure  followed  is  to  individually  sort  each  of  the  files 
and  then  merge  them  into  one  large,  ordered  file.  The  code  for  the  SORT  is: 

SORT,  VAR=DISC 
BYTES I ZE, 60 

FILE ,SORT=TAPEIN,OUTPUT=TAPEOUT 
FIELD, ROW (1 . 1,1, INTEGER) 

KEY, ROW (A) 

END 

The  merge  is  accomplished  in  a similar  fashion. 

With  an  angular  increment  of  1/2°  between  radial  scan  lines,  the 
description  of  the  radar  image  requires  1,440,000  points  (records).  The 
sort-in9  of  this  amount  of  information  is  a time-consuming  process  as  is 
evidenced  by  a SORT  running  time  of  1745  cp  seconds  (about  2 hours  wall 
clock  time)  and  a run  cost  of  around  347  dollars. 

10.  Program  RSS10 


10.1 

RSS10  uses  the  sorted  data  from  SORT  to  format  a plot  tape  for 

the  DICOMED. 

10.2 


The  DICOMED  screen  consists  of  a 2048  X 2048  grid  with  the 
intensity  of  each  grid  element  being  specified  by  a 6-bit  color  code.  In 
the  mode  of  operation  presently  being  employed,  plotting  is  done  sequentially 
in  horizontal  rows  beginning  at  the  top  of  the  screen.  Therefore,  in  order 
to  describe  a picture  it  ic  necessary  to  format  2048  records  of  205  words 
each,  e.g.,  one  record  for  each  row  consisting  of  2048  6-bit  color  codes 
packed  10  per  word. 

The  present  image  format  calls  fc*  the  radar  scene  to  be  displayed  as 
a circle  of  radius  1000  pixels.  Therefore,  the  first  step  is  to  generate 
11  blank  (all  white)  records  to  describe  the  top  margin  of  the  picture. 

The  actual  image  begins  in  row  12. 

For  each  of  the  2000  records  containing  image  information,  those  pixels 
actually  lying  in  the  radar-scene  circle  are  colored  black.  As  should  be 
evident,  not  all  points  within  the  circle  will  contain  image  information  - 
only  the  points  on  the  radial  scan  lines  contain  data  and  if  an  angular  spacing 
of  1/2°  is  used,  these  comprise  only  about  29%  of  the  total  number  of  possible 
image  points.  Therefore,  the  choice  of  the  "fill"  color  is  important  and 
black  is  chosen  since  it  will  not  introduce  any  correlation  error  when  the 
generated  image  is  compared  to  the  more  dense  output  of  a real  radar . 


RHMWnM 


HMH 

0 
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The  information  from  the  input  file  is  then  written  over  the  fill  color. 
Since  the  information  is  sorted,  taking  the  points  in  the  order  they  appear 
permits  the  records  to  be  filled  sequentially  from  left  to  right,  with  a change 
in  Y signalling  the  end  of  a given  line.  When  the  radar  circle  is  completed, 
a bottom  margin  of  11  blank  records  is  written  to  the  plot  tape. 

The  program  also  places  four  fiducial  marks  located  at  each  edge  of  the 
image.  A small  cross-mark  is  also  placed  at  the  center  of  the  image  to  mark 
the  target  location. 

10. 3 Output 

RSS10  outputs  the  DICOMED  plot  tape  which  marks  the  end  of 
the  simulation  process.  The  format  of  the  image  generated  by  this  tape  is 
presented  in  Figure  II.  11. 

10.4  Timing  and  Cost 

RSS10  requires  approximately  506  cp  seconds  execution  time 
and  about  20  minutes  wall  clock  time.  Run  cost  may  be  expected  to  be  in  the 
neighborhood  of  fifty- four  dollars. 
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III.  OPERATING  INSTRUCTIONS 


It  is  the  purpose  of  this  chapter  to  present  the  deck  structures 
required  to  operate  the  DRLMS  on  ETL's  CDC  6400  under  SCOPE  3.4.  As  is 
evident  from  the  preceding  discussions,  this  system  consists  of  a large 
number  of  independent  programs.  This  configuration  is  due  to  the  fact 
that  the  DRLMS  was  originally  written  for  a minicomputer  and  so  its 
structure  reflects  the  size  limitations  of  the  associated  hardware.  In 
order  to  simplify  and  speed  up  processing,  it  is  recommended  that  the 
system  be  run  as  two  dependent  strings  of  programs.  The  first  string 
consists  of  programs  RLMS1  through  RLMS6  and  may  be  viewed  as  a data  base 
preparation  subsystem.  Programs  RLMS7  through  RLMS10  may  then  be  run  as 
a dependent  string  to  create  the  radar  scene  for  the  desired  altitude. 

We  first  present  the  deck  structures  for  the  dependent  string  RLMS1 
through  RLMS6.  In  all  that  follows,  the  permanent  file  organization  is 
such  that  no  file  is  ever  purged.  The  user  may  wish  to  conserve  space  by 
altering  this  procedure.  Unit  numbers  must  be  assigned  as  indicated  al- 
though permanent  file  names  and  tape  VSN's  may  be  chosen  as  desired.  We 
emphasize  that  the  method  to  be  presented  is  only  one  possibility  and 
that  disc  storage  limitations  may  necessitate  the  use  of  more  space- 
effective  procedures. 

RSS1 : Here,  tape  7 is  the  Gerber  plot  tape  containing  the  plani- 

metry data. 

ETRS 1 , T200 , NTl , DABOO . 

TASK (TNET*****,PW*****,TRTS) NAME 
FTN. 

MOUNT (VSN=PK0007 , SN=BZRADAR) 

REQUEST ,TAPE8 , *PF,SN=BZRADAR. 

REQUEST , TAPE 7, NT , H D , S , US , NORING , VSN=EOO*  * * . 

FILE (TAPE7,MRL=1960, MB L= 1980, BFS=200) 

REWI N D ( TAPE 7 , TAPE8 ) 

LDSET, FILES=TAPE7 . 

LGO. 

CATALOG ( TAPE 8 , PFNl , SN=BZRADAR , I D=ET** *** , CY=0 1 , AC=RLMS 1 , RP=999 
TRANSF ( ETRS  2 , ETRS  3 , Ei RS  5 , ETRS6 , ETRS  7 , ETRS8) 

789 

FORTRAN  DECK 
789 

DATA  CARD  FOR  X,Y  OFFSETS 
6789 

R5S2: 


ETRS 2, CM200000 , T3000 , DAB01 . 

TASK (TNET*****,PW*****,TRTS) NAME 
FTN. 

MOUNT (VSN=PK0007,SN=BZRADAR) 

ATTACH ( TAPE8 , PFNl , SN=BZRADAR, I D=ET*  * * * * , CY=0 1 ) 
REQUEST, TAPE9, *PF , SN=BZRADAR 
REQUEST, TAPE20, *PF,SN=BZRADAR 
LGO. 
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CATALOG (TAPE9 , PFN3, SN=BZRADAR, I D-ET*** ** ,CY=02 , AC-RLMS2 ,RP=999) 
CATALOG (TAPE20 , PFN20, SN-BZRADAR, ID=ET*  ***  * ,CY=1 . AC=20UT, RP=999) 
TRANSF ( ETRS 3 , ETRS 5 , ETRS6 , ETRS 7 , ETRS 8) 

789 

FORTRAN  DECK 
789 

DATA  CARD  FOR  MAP  SCALE  PLANIMETRY  RESOLUTION  AND  REGION  SIZE 
6789 

RSS3:  RECALL  THAT  RSS4  IS  THE  SORT  ROUTINE  IN  RSS3. 

ETRS  3 , T1000 , DAB 02 . 

.ASK (TNET*****,PW*****,TRTS) NAME 
FTN. 

MOUNT (VSN-PK0007 ,SN=BZRADAR) 

ATTACH ( TAPE7 , PFN  3 , SN-BZRADAR , I D=ET  * * * * * , CY=02 ) 

ATTACH (TAPE20, PFN2 ,SN=BZRADAR, ID=ET***** ,CY=01) 

REQUEST, TAPE2, *PF,SN=BZRADAR 
LGO. 

FILE , TAPE1 , BT=C, RT=Z , FL=89 

FI LE , TAPE 2 , BT=C , RT=Z , FL=89 , CM- YES 

SORTMRG. 

CATALOG (TAPE2 , PFN4 , ID=ET*****, CY=03 , AC-RLMS3 , RP=999) 

TRANS  F ( ETRS  5 , ETRS6 , ETRS  7 , ETRS  8 ) 

789 

FORTRAN  DECK  FOR  RLMS  3 
789 

SORT/MERGE  DIRECTIVES 
6789 

RSS5: 

ETRS 5, Tl 000, DAB 3. 

TASK (TNET***** ,PW***** ,TRTS) NAME 
FTN. 

MOUNT (VSN-PK0007 , SN=BZRADAR) 

ATTACH ( TAPE  7,PFN4,SN=BZRADAR,ID=ET*****,CY=03) 

ATTACH (TAPE20 , PFN2 , SN=BZRADAR , I D=ET*  * *** , CY=01 ) 

FI LE , TAPE 7 , BT=C , RT=Z , r L=89 . 

REQUEST, TAPE10, *PF,SN=BZRADAR. 

LDSET ( FI LES=TAPE 7 ) 

LGO. 

CATALOG (TAPE 10 , PFN6 , SN-BZRADAR, ID=ET** ** * , CY=05 , AC=RLMS5 , RP=999) 
TRANSF (ETSM6 , ETRS7 , ETRS8) 

789 

FORTRAN  DECK 
6789 
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RSS6:  Here,  TAPES  is  the  UNAMACE  elevation  tape. 

ETRS6 , T2000 , CM2  30000 , NTl , DAR04 . 

TASK (TNET*****,PW*****,TRTS) NAME 
FTN. 

MOUNT ( VSN=PK0007 , SN=BZRADAR) 

ATTACH  (TAPE20 , PFN2  , SN=BZRADAR,  I D=ET***  **  , C'Y  = 01 ) 

REQUEST  (TAPE8,  NT, flD,  S , NORING=VSN=EOO**  *) 

REQUEST (TAPE9, *PF,SN=BZ RADAR) 

LGO. 

CATALOG (TAPE9,PFN7,SN---BZRADAR,  I D=  ET*** **. CY=04  , AOPROG6  . RP=999) 

789 

FORTRAN  DECK 
6789 

This  completes  the  data  base  preparation  stage  of  the  DRLMS.  We 
emphasize  again  that  the  preceeding  programs  need  be  run  only  once  for 
each  target  area.  We  now  proceed  to  the  image  preparation  stage  con- 
sisting of  programs  RSS7  through  RSS10. 

RSS7: 

ETRS7 ,T700 , DAB05 

TASK (TNET*****,PW*****,TRTS) NAME 
FTN. 

MOUNT (VSN=PK0007 ,SN=BZRADAR) 

REQUEST, TAPF.6,  *PF,SN=BZRADAR. 

REQUEST , TAPE  12 , *PF,SN=8ZRADAR. 

ATTACH ( TAPE  3 , PFN6 , SN=BZRADAR , I D=ET*  * * * * , CY=05 ) 

ATTACH (TAPE! , PFN7 ,SN=BZRADAR , ET* **** , CY=04) 

ATTACH (TAPE20,PFN2.SN=BZRADAR,ET*****,CY=01) 

LGO. 

CATALOG ( TAPE6 , PFN8 , SN=BZRADAR, 1 D=ET*  * * * * , CY=04 , AO RLMS 7 , RP=999 ) 
CATALOG (TAPE1 2, PFN9 , SN=BZRADAR, ID=ET***** ,CY=04 , AC=MAPDATA , RP=999) 
TRANSF (ETRS8) 

789 

FORTRAN  DECK 
789 

$PARAMS 
DATA  CARD 
$END 
6789 

RSS8 : 

ETRS8 , T4500 , DAB06 

TASK (TNET*****,PW*****,TRTS) NAME 

FTN. 

MOUNT ( VSN=PK0007 ,SN=BZRADAR) 

ATTACH (TAPE6,PFN8,SN=BZRADAR, I D=ET*****,CY=04) 

ATTACH (TAPE12,PFN9,SN=BZRADAR, ID=ET*****,CY=04) 

ATTACH (TAPE 20 , PFN2 , SN=BZRADAR, ET***** , CY=01 ) 

REQUEST , TAPE4 , *PF , SN=BZRADAR. 

LGO. 

CATALOG (TAPE4 , PFN10 , SN=BZRADAR , ID=ET***** , CY=32 , AOPROG8 , RP=999) 
789 


RSS9 : 


The  following  must  be  rut  once  for  each  desired  radar  altitude. 
Note  that  the  output  of  1SS9  is  placed  on  a private  pack.  This 
is  duo  to  space  limitatitns  on  ETL's  system  disc. 


11 

fi 


li 

ll 

u 


n 


r ] 

1 1 


i:TRS9 . T 3500 , DAB00 

TASK (TNET***** , PW***»* .TRTS) NAME 

FTN. 

MOb NT (VSN=PK0G07 , SN=BZRADAR) 

ATTACH i TAPE6 , PFN6 , SN=BZRADAR, ID=ET*****, CY=04 ) 

ATTACH (TAPE4 , PIN10 ,SN=BZRADAR, ID=  ET*»*** ,CY=32) 

ATTACH ( TAPE2C , PFN2 , SN=BZRADAR, I D=  ET*  * * * * , CY=01 ) 

REQUEST (TAPE  3 . *PF , SN=RZRADAR) 

REQUEST ( TAPE 10 , *PF, SN=BZRADAR) 

REQUEST ( TAPE 1 2 . *PF . SN=BZRADAR) 

REQUEST (TAPE14 , *PF, SN=BZRADAR) 

FILE (TAPE  3, BT-K , RB=  31 , RT=F , FL-20 , MBL=620) 

FILE (TAPEl 0 , BT=K , RB=  31 , RT=F , FL=2U, MBL=620) 

FILE (TAPE 12 , BT=K , RB=  31 , RT-F, FL-20 , MBL=620) 

FILE (TAPE14 , BT=K , RB=31 , RT=F, FL=2i  >,  MBL=620) 

LDSET (FITES -TAPE 3/TAFE1 0/TAPEl 2/TAPE 1 4 ) 

LGO. 

CATALOG ( TAPE  3 , PFNSN , I D=ET*  **** , SN=BZRADAR . CY=01 , AC=9O0T , RP=999) 

CATALOG (TAPE 10, PFNSN, 1 D=ET***** , SN-BZRADAR, CY=02 , AC=90UT, RP=999) 
CATALOG  t TAPE 1 2 , PFNSN , I D=ET*  * * **  . SN=BZRADAR, CY=0  3 , AC=90UT , RP=999) 
CATALOG (TAPF14 .PFNSN, ID=ET* * * ** , SN=BZRADAR, CY=04 , AC=90UT, RP=999) 

TRANSF (ETSORT . ETRS10) 

789 

FORTRAN  DECK 
6789 

ETSORT:  Due  to  space  limitations,  the  SYSTEM/SORT  scratch  files  are 

assiqned  to  the  private  pack.  These  are  denoted  by  the  ZZZZZ 
prefix. 

ETSORT , T5000 , LAB 01 . 

TASK  ( TNF.T*  * * * * , PV7*  * * * * , TRTS ) NAME 
MOUNT ( PK0007 , SN=BZRADAR) 

ATTACH (TAPE  3, PFNSN , ID-ET**  *** , SN-BZRADAR , CY=01 ) 

ATTACH  (TAPF.10, PFNSN,  ID=ET*****,SN=BZRADAR,CY=02) 

ATTACH (TAPE 12, PFNSN, ID=ET***** ,SN=BZRAUAR , CY=03) 

ATTACH ( TAPEl 4 , PFNSN , I D=ET  * * * * * , SN=BZRADAR , CY=04 ) 

REQUEST (ZZZZZ1 A, SN=BZRADAR) 

REQUEST  ( ZZZZZ  1 B , SN=BZRADAR) 

REQUEST ( ZZZZZ1 C . SN=BZRADAR) 

REQUEST! ZZZZZ i D,SN=BZRADAR) 

REQUEST ( ZZZZZ1E ,SN~BZRAPAR) 

REQUEST ( ZZZZZ1F , SN=BZRADAR) 

REQUES  T (TAPE 2 ,SN=BZRADAR) 

REQUEST (TAPE 11 ,SN=BZRADAR) 

REQUEST (TAPE  1 3 , SN=BZRADAR) 

REQ7JEST  (TAFE15  ,SN=BZRADAR) 

REQUEST (TAPE27, *PF, SN-BZRADAR) 

FILE (TAPE3 ,BT-K, RB=  31 , RT=F , FL=20 , MBL=620) 

FILE  (TAPE2,BT  K,RB=31  ,RT=F,FL=20,MBL=620) 

SORTMRG. 
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( 

1 


i. 


F J LE  lT/‘l  K10,BTK,RB  3]  . RT=F , FL=20 , MBL=620) 
FI li  (TAi  II, PT~K , RB=  31 , RT- F , FL -20  MBL=620) 

SORT* IK  > . 

FI  LE  ( l’Af  F 1 2 , UT-K  , RB~  31 , RT=F, FL-  20 , MBL=620) 
FI1.E(T,"'  :l  3.BT  K , RB - 31  . RT^F , FL=20 , MBL~620) 

SOr’TMK. . . 

FI  I K (TAI  111 4 , RT-  K , RB=  31 , RT-F , FL=2C , MBL=620) 
FI  1.E  (TAPF.l  5 , BT  - K , RB=  31 , RT=F , FL=20 , MBL=620) 

sor^mr;;. 

FILE  ( i A1'1.27,BT=K.RB=31  , RT=F , FL=20 , MBL=620) 
SORTMRC  . 

CA  TAU  v ; I TAP.  2 7 , F I MALE , SN=  BZRADAR , I D=ET  * * * * 
TRANSF  li-.TRSlO) 

"'89 

SORT  DIRECTIVES  FOR  FIRST  SORT 
789 

SORT  DIRECTIVES  FOR  SECOND  SORT 
7h9 

SORT  DIRECTIVES  FOR  THIRD  SORT 
7P) 

SORT  DIRECTIVES  FOR  FOURTH  SORT 
789 

MERGE  DIRECTIVES 
6789 


, AC=PLOTPE . RP=999) 


It  should  be  clear  that  all  four  SORT  routines  are  the  same  except  for  the 
file  names. 

3SS 10: 

ETRS10 . Il 500,TFl,DAB02. 

TASK (TNET*****,PW*****,TRTS) NAME 
FTN. 

MOUNT ( PK000? , SN=BZRADAR) 

ATTACH ( TAP  E 1 . FI  NA LE , I D=ET  * * * * * , S N=  BZRADA R ) 

FUJI  (TAPEl , B T=K , RB=  31 , RT=F  , FL=20 , MBL=620) 

REQUEST, TAPE  3 , L, RING , VSN=E00*** . 

LDSLT (I  I LES=TAPE1 ) 

LGO. 

789 

FORTRAN  DECK 
6789 


The  output  from  RSS10  is  the  plot  tape  for  the  DICOMED. 
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FIGURE  IV. 1 - PROGRAM  RSSi  FLOWCHART 
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